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; Abstract 

The goal of this paper is to develop the formalism of the two-particle irreducible (2PI) 1 (or Cornwall 
- Jackiw - Tomboulis (CJT) effective action (EA) in a way appropiate to its application to non 



00 



o 

Oh 



- 1—1 

X 



equilibrium gauge theories. We hope this review article will stimulate new work into this field. 

1 Introduction 



The goal of this paper is to develop the formalism of the two-particle irreducible (2PI) [T] (or Cornwall - 
Jackiw - Tomboulis (CJT) [2]) effective action (EA) in a way appropiate to its application to non equilibrium 
gauge theories. The usual formulation of the 2P1EA cannot be extended to these theories because of the 
special features of gauge invariance. There are two difficulties in particular which require consideration, 
namely, the existence of constraints linking the Schwinger functions of the theory among themselves [3], and 
the peculiarities of the gauge - fixing procedure. The former are expressed by the so-called Takahashi - Ward 
or Slavnov - Taylor identities, which in turn derive from the Zinn - Justin equation (see below) 0]. The 
latter raises the issue, exclusive to gauge theories, of the gauge fixing dependence of theoretical constructs 
0. A clear understanding of this issue is essential to the physical interpretation of the theory. 
Oh, Non equilibrium quantum field theory [5] has evolved in the last few years mostly in the context of 

non - gauge field theories. In these studies, two tools have proved extremely valuable, namely the closed 
- time - path (CTP), IN-IN or Schwinger - Keldish formalism Q, and the 2PI or CJT EA. The former 
allows us to study the causal evolution of quantum fields, in contradistinction to the IN -OUT or mixed 
advanced and retarded boundary conditions appropiate to the study of scattering processes. The latter 
provides a comprehensive framework where various non - perturbative approaches may be most efficiently 
implemented. The need to go beyond simple - minded perturbation theory is an universal feature of high 
energy nonequilibrium processes in non linear theories. 

When we survey the literature on gauge theories, we find that the one - particle ireducible (1PI) IN - OUT 
EA is a well developed tool which has found its way into most modern textbooks The CTP formulation 
of gauge theories, although less widespread, has also been the subject of several important investigations 
and may be considered well understood |1(J1 lllj . The 2PIEA, on the other hand, has only recently come 
under study E| ■ We hope this review article will stimulate new work into this field. 

Of course, the subject of non equilibrium gauge theories is so vast that it becomes impossible to make 
progress without some well defined self - impossed limitations from the outset. We will restrict ourselves to 
Yang - Mills and to non linear abelian theories such as QED and SQED. We shall make no explicit attempt 
to discuss gravity, form fields or string theories |14|. 

These self - imposed limitations in aims are correlated with some necessary a priori technical choices. 
We shall discuss only the path integral Fadeev - Popov quantization of gauge theories. Although we shall 
use BRST invariance at several stages, we shall not apply methods such as BRST or Batalin - Vilkovisky 
quantization, which really come on their own only in more demanding applications [5]. We shall use DeWitt's 
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notation and are deeply in debt to DeWitt's insights ^5], but we shall not use the gauge - independent 
formulation of DcWitt and Vilkovisky J^](on this subject, see the discussion in 17), nor more recent 
developments by DeWitt and collaborators |18| . 

When gauge symmetries are unbroken, there are no preferred directions in gauge space, and all back- 
ground fields will vanish identically. Therefore, the only degrees of freedom in the 2PI formalism shall be the 
propagators or two - point functions. Also, there will be no need to distinguish between the usual and the 
DeWitt- Abbott gauge invariant EA nor to introduce gauge fixing conditions appropiate to the study 
of broken gauge theories, such as the i?£ family of gauges [H]- We shall only assume that the gauge fixing 
condition is linear on the quantum fields. On the other hand, we shall be completely general regarding group 
structure, matter content, (linear) gauge fixing condition and gauge fixing parameter. 

The main results of this paper are formulae l|78l) and l|84|) providing the definition of the 2PIEA for gauge 
theories. We shall then use this construction to discuss the Zinn-Justin identities and the gauge dependence 
of the propagators. These are in some sense known results, and we include them to help the reader to connect 
the 2PI to earlier formulations of non equilibrium field theories, and to better appreciate its power. 

The paper is organized as follows. Section 2 is a review of the path integral approach to gauge theories, 
including the IN-OUT effective action and the Kugo canonical formulation [201 • This section stablishes 
our notation, and sets the standard for the new developments that follow. Section 3 contains the main 
results, including the construction of the 2PIEA and the proof of its structure as the sum of 2PI Feynman 
graphs. In Section 4, we use the 2PIEA as starting point for the discussion of the Zinn - Justin identity. In 
Sections 5 and 6 we use the ZJ identiry to study the gauge dependence and the structure of the propagators, 
respectively. 

We have gathered in the Appendix some relevant formulae concerning Grassmann calculus. For more 
details, we refer the reader to the monographs by Berezin j2J, DeWitt [22] and Negele and Orland [22| • 

2 Path Integral approach 

2.1 Gauge theories 

A gauge theory contains "matter" fields ip such that there are local (unitary) transformations g which are 
symmetries of the theory. The g's form a non abelian (simple) group. Infinitesimal transformations may 
be written as g = exp [is] , where the hermitian matrix e may be expanded as a linear combination of 
"generators" e = £ A Ta- The generators form a closed algebra under commutation 



The structure constants C AB are antisymmetric on A, B and satisfy the Jacobi identity. 

Gauge invariance of kinetic terms within the Lagrangian means that derivatives are written in terms of 
the gauge covariant derivative operator = — iA^. The conexion A M = A^aT a transforms upon an 
infinitesimal gauge transformation as 



[T A ,T B ]=iC c AB T c 



(1) 




A^ + D M e 



(2) 



where 



D^e = 9 M e - % [A M ,s] 

Covariant derivatives do not commute, but their commutator contains no derivatives 



(3) 



[D^D V ] = -IF) 



(4) 



where the field tensor 



dfj,Av di/A^ — i [A^, Ajy] 



(5) 



Upon a gauge transformation 



F, 



F„ v + i [e, Fpu] 



(6) 
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therefore the object 



Ag 2 



(7) 



is gauge invariant. This is the classical Lagrangian density for the gauge fields, g being the coupling constant. 
The total action S = S + S m , where 



So = J d d x (_T)TrW^ 
and S m is the gauge invariant action for the matter fields. 



(8) 



2.2 DeWitt's notation 

We may drop the distinction between gauge and matter fields, and consider a theory described by a string 
of fields <b a invariant under infinitesimal transformations 



j a rria r n _A 



Time 1 



(9) 



The commutation rules are the statement that the commutator of two gauge transforms is also a gauge 
transform, namely 



The classical equations of motion read 



8S 



= 



and because of gauge invariance we must have the identity 



SS 



(10) 



(ii) 



(12) 



2.3 The vacuum to vacuum amplitude 

In the quantum theory, we expect the vacuum to vacuum amplitude to be given by the IN-OUT path integral 



D<t> e 



iS 



(13) 



However this integral counts each history many times, and is ill defined. 

To cure this problem, let f A be functionals in history space which are not gauge invariant. This 
means that, given a history <p a such that f A [<f) a ] — 0, and an infinitesimal gauge transform such that 
f A [<fi a + 5<p a ] — too, the gauge transform must be trivial; in other words 



Dot 



Sf A 



^0 



(14) 



Now let us call </> [e] the result of applying a gauge transform parameterized by e to the field configuration 
b. Then we have the identity 



/ 



Be Dct 



S[f A [<f>{e}}~C A ] =1 



where C may be anything, and we can write 



Z = De Dcf) Det 



[e]]TS[4>[e]] 



S [f A [ft [s]] - C A ] e lS ^ 



(15) 



(16) 
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Of course, S [0] = 5 [(f) [e]] ,and 



D<f> [e] = D(j> { 1 



so, provided 



we find, up to a constant 



Tr- 



= o 



Z= D<f> Det 



c 



Al e iSf[, 



(17) 



(18) 



(19) 



Since the C A are arbitrary, any average over different choices will do too. For example, given a suitable 
metric we may take the Gaussian average 



DC A e -( l / 2 t)c A c A 

Integrating over £ and after a Fourier transform we find 

~5f A 



Z 



D(j)Dh A Det 



mm 



exp < i 



S [0] + h A f A [</>] 



i ' ^-h A h A 



(20) 



(21) 



h A is the Nakanishi - Lautrup (N-L) field |2U, an d £, the gauge fixing parameter. 
2.4 Ghosts 

We may write the determinant as a functional integral 



Z = Duj b D X AD(j)Dh A exp <^ i 



S [</,] + h A f A [</»] 



h A h A + i XA A A 



(22) 



(23) 



The uj b , xa are independent c- number Grassmann variables, namely the ghost and antighost fields, 
respectively. Following Kugo and Ojima and unlike Weinberg [S], we have included a factor of i in the 
ghost Lagrangian, which is consistent with taking the ghosts as formally "Hermitian" and demanding the 
action to be "real" . 

We assign "ghost number" 1 to uo B , and —1 to \A- 

2.5 BRST invariance 

We may regard the functional 



S eff = S[ct>]+h A f J 



Ik* 
2 



h A h< 



tXAA J 



(24) 



as the action of a new theory, built from the original by adding the N-L, ghost and antighost fields. By 
construction, this action is not gauge invariant in the original sense. However, let us consider a gauge 
transform parameterized by 0lo b , where is an anticommuting "constant", namely 



5(f> a = 9T% [cj)} uj a 
Observe that, keeping the other fields invariant for the time being 

Sf A [0] = 6A A 



(25) 
(26) 
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5A A = QJ% (i T% [4>] [4>] lu b lu c + 9f%T% [4>\ >/3 Tq [0] lu c u; b (27) 
Since the ghosts are Grassmann, this becomes 

6A A = l -Qf A a T% [</>] C d bc u c uj b (28) 
These results suggest extending the definition of the transformation to 

Sh A = (29) 
S X A = iQh A (30) 

Suj d = hc D CB uj c uj B (31) 
Then S e ff is invariant under this"BRST" transformation. Let us define the operator f2 

n [X] = ±5X (32) 
The operator f2 increases "ghost number" by one. Obviously f2 2 [Ha] = ^ 2 [xa] — 0, but also 

n 2 [4> a ] = T% [</>} >/3 T B [0] to B cu A + \T% [<j>] C a cb u c lo b = (33) 

r>2 [ w ^] = [C c £f /// + C b ef lu c lu e lo f ] = (34) 

These properties imply that actually fi 2 = tout court (see |S]). Also, observe that 

S eff = S + Q[F] (35) 



where 



is BRST invariant, and 



S = S[<p] (36) 



^ = -*XA + (37) 

(recall that [F] = -i (17 [xa] {/ A [0] + \£,h A ) - x^ft [/ A [<£[])). Also, observe that, provided 

C\ B = (38) 
the functional volume clement is also BRST invariant. 

2.6 Gauge invariance - gauge independence of the vacuum to vacuum amplitude 

It follows from the above that any gauge fixing dependence (that is, dependence on the choice of the gauge 
fixing condition f A , gauge fixing parameter £ or the metric used to raise indexes in the N-L field) may only 
come from a dependence upon changes in the functional F. Any such change induces a perturbation 

5Z = i J Duj b 'D XA D(f>DhA & [SF] exp {iS eff } (39) 
Now, call X r the different fields in the theory. Then 

n [SF] = (-i)^ +1 SF. r n [X r ] (40) 
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where g r is the corresponding ghost number. Integrating by parts (see )26|). and provided the surface term 
vanishes, we get 



SZ = i J Duj B D XA D<t>Dh A 6F l-^exp{iS e ff}n[X r ] 



(41) 



But the bracketts vanish, because of BRST invariance of S e ff and because ft [X r ] is divergence - free. 
Therefore the physicality condition is that the flux of any vector pointing in the direction of ft [X r ] over the 
boundary of the space of field configurations must vanish. 

This shows by the way that F could be any expression of ghost number —1, since S e ff must have ghost 
number zero. 

2.7 In-out, vacuum effective action 

Let us write the generating functional for connected Feynman graphs 

exp{iW} = J Du B D XA D<j)Dh A exp{i [S eff + J r X r }} (42) 
and define the in-out effective action as its Legendre transform 

T [X a ] = W[J\- J r X r (43) 

Performing a change of variables within the path integral corresponding to a BRST transformation, the 
generating functional cannot change, and therefore 

J Du B D XA D(j)Dh A exp {i [S eff + J r X r }} J r Q [X r ] = (44) 

but 

J r = -I\ r (45) 

and so we obtain the Zinn - Justin equation 

r r (fi [X r ]) = (46) 

where 

(R) = cxp{-iW} J Duj b D X AD(j)Dh A exp{z [S eff + J r X r }} R (47) 

On the other hand, if we simply change the functional F in S e ff by an amount SF, we get, holding the 
sources constant 

SW = exp {-iW} J Du B D XA Dcj)Dh A VL [8F] exp {i [S eff + J r X r }} (48) 
and repeating the previous argument 

§W = iT ir (SF Q [X r ]} (49) 

Of course, if we hold the background fields constant, the sources will change by an amount 5J r . However, 
this extra variation does not contribute to the Legendre transform, and so 

ST = iT >r (SF n [X r ]) (50) 



G 



2.8 The canonical approach and non - vacuum states 

In this section, we shall consider the concrete case where <p a = A A , f A = S B d^ and Tg = Sgd^ + C A B A^. 
We can write S e ff explicitly 



S eff = J d 4 x l^F A ^F A ^ - d^h A A» A + ih A h A - id^XA [&id» + iCg B A%\ w s | 



(51) 



If we take A Aa (a — 1,2,3), h A , xa and uj a as canonical variables, then we may identify the corresponding 
momenta (201 

p Aa = ]_ F A a0 . p A = _ A A0, p A = _ % ^ + /; p uA = id oX A (52) 

and impose the ETCCRs 

\px r ,X% = -i5 r s (53) 

where we use anticommutators for ghost fields and momenta, and commutators for all other cases. 
The BRST invariance of S e f / implies the conservation of the Noether current 

r- a m%& (54 » 



We define the BRST charge 



n = J d 3 xn[x r ] PXr (55) 

This is the generator of BRST transforms, since 

sx r = en [x r ] = i [en, x r ] (56) 

(since 9 is Grassman, we use commutators throughout). Then n 2 = 0. 
S e f f is also invariant upon the scale transformation 

c^-eV, X B ^e-\ B (57) 

The corresponding generator 

Q = J d 3 x {u B p uB - XAP A } (58) 
is the ghost charge. Ghost charge is bosonic, so [Q, Q] = 0. On the other hand, n has ghost charge 1, so 

i[Q,Cl]=n (59) 

Both Q and en commute with the effective Hamiltonian. 

Observables are BRST invariant, and so they commute with 6T2. Physical states are also BRST invariant, 
therefore annihilated by fl. Physical states differing by a BRST transform are physically indistinguishable, 
in the sense that they lead to the same matrix elements for all observables, and so we may add the condition 
that a physical state \a) is BRST-closed (f2 \a) = 0) but not exact (there is no \0) such that |a) = ft \(3)). 

One - particle unphysical states come in tetrads |20]. Let \N) be a one - particle unphysical state with 
ghost charge N, i. e., iQ \N) — N \N) (any further BRST invariant quantum numbers are irrelevant to the 
argument, and shall be omitted). Call \N + 1) = O, \N) . Then fi \N + 1) = 0, and this implies that |JV + 1) 
has vanishing norm, since (N + 1 \N + 1) = (N + 1| Q, \N) = 0. Let \-N - 1) be the state "conjugated" to 
\N + 1), in the sense that (-N — 1 \N + 1) = 1 (since iQ is hermitian, this state must have ghost charge 
-N- 1). Now (-N — 1| fl \N) = 1, and so \-N) = fl \-N - 1) is conjugated to \N) . Observe that, without 
loss of generality, we may assume that N is even. 

Now let a/ M be the corresponding creation operators, namely \M) = a} M |0) (M — N, —N, JV + 1 or 



N — f). Then a^ N+1 = Cl, ajy and aJ_ N = |o, aLjv-i }■ ^ow ^ 2 = ^' so |^' a lv+i} 



t 

N 



0. 
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o—N, a_ N 



0; the 



Because |iV+l) and \—N) have vanishing norms, we must have ■fajv+i, ojv+i j 

conjugacy relations imply |a_^v_x, ajy + i j = a-Af 7 ajv = 1- The Oji^'s are destruction operators, but ajy 

destroys the |— iV) (rather than the \N)) state, oat+i destroys |— N — 1) , etc. 

It follows that the projection operator Pi over the subspace of states with exactly one unphysical particle 
may be written in terms of the projector P over physical states as 



But then 



Pi 



l — N 



Pa N + alfPa-N + a f _ N _ 1 Pa N+1 + a^ N+1 Pa_ 



N-l 



(60) 



Pi 



9 



9 



9 



''-N-l 



Pa 



N 



'-JV-1 



Pa 



N 



-N- 



-N-l 



Pa 



N 



-N- 



a' N Pa-N-i 



9PaN + a^Po- 



-N 



l -N-l 



Pa 



N+l 



a' N Pa-N-i 



1 P9aw + ajyPa- 



iV 



-JV- 



_ t Pa 



N+l 



I f 4 



ajyPOa-Ar-i 



_ 1 P9ajy + ajv-P a -N-i} — a -AT-i-P [^i a iv] ^ a] v P9 



|n, jV _ 1 Pa JV + ajyPa-jv-i j 



(61) 



Repeating identical arguments for the spaces with n unphysical particles, we conclude that the projector 
P' orthogonal to P has the form P' = {9, R}, where R is some operator with ghost number —1. 

We may now deal with the construction of statistical operators in gauge theories. In principle, a physical 
statistical operator should shield nonzero probabilities only for physical states, and so it should satisfy 
p = Pp = pP. This is a much stronger requirement than BRST invariance [9, p] = 0. So, given a BRST 
invariant density matrix p, we ought to define the physical expectation value of any (BRST invariant) 
observable C as 



(C) phys =Tr[PpC] (62) 

However, Kugo and Hata [2HI (KH) have shown that the same expectation values may be obtained by 
using the statistical operator e~™Q p . The key to the argument is that the commutation relation [iQ, 9] = 9 
implies that, if \N) is an eigenstate of iQ with eigenvalue N, then 9 \N) has eigenvalue N + l. It follows 
that {e~ 7rQ , 9} = 0, since e _7rQ = e™ {iQ \ We then find that, for any BRST invariant observable C 

(C) phys = Tr [PpC] = Tr [Pe^pC] = Tr [e^ Q pC] - Tr [{9, R] e^ Q P C] (63) 

We must show that the second term vanishes, and this follows from {e~ 7rl 3, fi} = and [9, pC] = 0. 

This suggest to define the expectation value (C) of any observable as (C) = Tr [e _7r< 2/oC] . Of course, this 
agrees with the physical expectation value only if C is BRST invariant. For example, the partition function 
computed from e _7r( 5p agrees with the partition function defined by tracing only over physical states, but 
the generating functionals obtained by adding sources coupled to non - BRST invariant operators will in 
general be different. 

The advantages of the Kugo - Hata ansatz are clearly seen by considering the form of the KMS theorem 
appropiate to the ghost propagator. Let us define 

G%(x,x') = (P[ X a A (x)u; B (x')]) (64) 
where P is the usual (CTP)-ordering operator. Then 

G a 1 b (x,x')=( X a(x)u Jb (x')) (65) 
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G a 2 b (x,x') = -(cj b (x , ) X a(x)) (66) 

(observe the sign change, associated to the anticommuting character of the ghost fields). The Jordan prop- 
agator is defined as G — G 21 — G 12 



If we omitted the K-H e factor, we would reason, given p - 



e 



-0H 



G 2 2 B (x, x 1 ) « Tr [e-0 H X A (x) to B (x')] = Tr [ XA (x + i/3) e^ H io B (x')] = -G\\ (x + i/3, x') (67) 

Therefore G 2 AB (oj) = —e^ u G A 2 B (to) , leading to a Fermi - Dirac form of the thermal propagators. 
This reasoning is incorrect. The proper way is 

G 2 2 B (x, x') = Tr [e-^XA (x + i/3) e~^ H u B (x')] = G 12 B (x + i/3, x') (68) 
So G 2 AB (uj) = e^GW (uj), which leads to the Bose - Einstein form. 

Let us observe that in the path integral representation, the K-H factor does not appear explicitly, but 
only changes the boundary conditions on ghost fields from anti-periodic to periodic. 

We conclude that in this formalism, unphysical degrees of freedom and ghosts get thermal corrections, 
both being of the Bose - Einstein form, in spite of the ghosts being fermions (for which reason ghost loops 
do get a minus sign). For an alternative formulation, see 

3 The 2PI formalism applied to gauge theories 

We can now begin with our real goal, namely, the application of the 2PI CTP formalism to gauge theories. 
We shall proceed with a fair amount of generality, only assuming that the gauge condition is linear, and also 
the gauge generators T% [0] = Tg A + T^ A0 ^ . 
The classical action is given by 

S eff =S[<f>] + h A f A [4>] + h A h A + ixaj^ [4>] T B [4>] u B (69) 
To this we add sources coupled to the individual degrees of freedom and also to their products 

X r J r + ^X r K rs X s = j a x a + 9 U X U + ^n ap x a x p + ^a uv 9 u 6 v + 9 u i> ua x a (70) 

where x a represents the bosonic degrees of freedom {<fi, h) and 9 the Grassmann ones (ui, x), and we 
introduce the definition K QU = — K UQ . Observe that j, k and a are normal, while A and ip are Grassmann. 
a is antisymmetric. 

We therefore define the generating functional 



e lW = I DX r exp <j * 



&eff + X r J r + —X 1 K. rs X 6 



(71) 



Note that the information about the initial state is implicit in the integration measure, and will reappear 
only as an initial condition on the equations of motion. 
We find 

W^t = X r (72) 

,J r 
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isrns 



W-— = ^— X r X s + G rs \ (73) 
oK rs 2 

where we introduce the bookkeeping devise 9 r = (— l) 9r , where q r is the ghost charge of the corresponding 
field, and 9 rs = (-1) 9 *" 33 . 

We define the Legendre transform 



9 



r = w - x r j r - -x r K rs x s 



-G rS K r 



whereby 
S 



sx r 



-j r - e rs x s 



sx r 

S 



+ > 



w 



SK S 



5X r ~ 



-K s 



-J. 



SX 

1 f) S ^f)^ 

■~e sr x s K sr - - — e rs e rt G st 

2 2 



- -K rs x s - -e rs x s 

2 2 



SX 



5X r 



-K, 



Qstat 



-X'X* 



SX 



Jr 2 ^-rs-X 



SX 



X' 



-J r — -K r cX s — -^ r X s K^ 
2 2 



Now observe that K, r = 6» r 6> s 6>'' s K r ,. In the end 



r — — t/ r — ic^o^x" 5 



In the same way 
6 -V = 



<5G r 







f 5 j1 












|_£G rs J 



_ O rt st x t 



Qtqgq 



_QTt QS t x t 



rtQStgrqQSq 
5 



w 



5K 



SG r 



K 



tg 



tg 



X 9 



(5G*- 



-K 



tg 



tgrqgsqQtq 



9 tq q - - 

= e rt 6 st e rq e sq x'x 9 

2 



-K, 



SG r 



-K 



tg 



<5G 
1 



-K 



tg 



2QrtQst X t 

2 



<SG' 



-K 



6g 



3.1 The 2PIEA 

In order to evaluate the 2PIEA, let us make the ansatz 



r = S [X r ] + -0 sr 6 s G rs S rs - - lnsdet [G rs ] +T 2 - -6 s G rs G R l 



B.sr' 



where 



SX r 



-s 



SX S 



and S is the classical action Eq. (|69|l . 

The corresponding ansatze for the sources are 



s 



■r 



-S r + Jl 



K, 



-26» sr (9 s 



SG r 



SX r 

= — S rs + *G ira + K2,- s 
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The generating functional 



W = r + X r J r + -X r K rs X s + -^G rs K rs 

= & [X r ] + \e sr 9 s G rs S rs - l - lnsdet [G rs ] + T 2 - ^G^G^ + X r [-S. r + J 2r ] 

1 - sr s 

--X r [-S„ + iGll s + K 2rs ] X s + —G rs [S r , + iGll s + K 2rs ] 

= S [X r ] - l - lnsdet [G rs ] + T 2 -X r (s, r - J 2r ~ \ (S rs - iG^ s - K 2rs ) X' 

nsr as 

H ^ — 2rs 



(82) 



so 



(sdct [G rs ]) 1/2 exp [«r 2 ] 
= J DX r cxp i [S (X) + X r (-5 r + J 2r ) 

+ (-S„ + iG~ L l s + K 2rs ) X s - X r (-S„ + iG£* a + K 2rs ) X s 

- ( - 1 - i i 1 
—5 [X r ] + X r f S r — -S rs X s + 2^i,rs^ S — ^ 2r + 2"^" 2r;s ^' 



-G rs r r 6» s K 2rs 



(83) 



or else 



jr 2 



where 



(sdet[G rs ])" 1/2 / DSX r expi 



AS + -8X r G- L l s 5X s + 5X r J 2r + [SX r SX s - G rs ] K 2rs 



(84) 



AS = S [X r + SX r ] - S [X r ] - 5X r S. r - -5X r S rs 5X s 



(85) 



T 2 is the sum of 2PI vacuum bubbles in a theory with free action iG^ rs and interacting terms coming from 
the cubic and quartic terms in the development of S around the mean fields. 

In spite of appearances, the new term 9 s G rs G^ r is a constant. To see this, parametrize 



G r 



leading to 



H a P N' 
N U P M 



VOL 

uv 



(86) 



r- 1 



H a p N va 
N u0 M uv 



(87) 



Then, because they are inverses, we must have 

H at3 H fil + N va N vl = 5° 



H a ^N vP + N ua M uv = 



N u0 Hr^ + M uv N vy = 
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N u0 N vli + M UW M WV = 6% 
and therefore 

s G"G^ r - {H af} H [ja - N va N va ) + (N u ?N u p - M UV M VU ) 

= 5Z~2N va N va -(5l-2N uP N u p)=5Z-5 u u (89) 

independent of G rs . It may therefore be discarded. 

3.2 The reduced 2PIEA 

Let us now investigate the Schwinger-Dyson equations 

5 -r = 



sx 

8 



5G r 

From Eq. (J7HJ) we get 



■r = (90) 



s -^s [x r ] + 1^w^g^^_s m + ^r 2 = o 



6_ 

'SG 

The second set of equations may be rewritten as 



9 sr e s s rs - iew* (G^) rs + 2— r 2 = o (91) 



6 

'SG 

and finally as 



S rs -z9 r (Gl 1 ) s +2—T 2 = (92) 



G 

The classical action is given by Eq. I|tj9fl . If we expand X r = X r + SX r , then the quadratic terms are 



S rs -i(Gi 1 ) rs + 20"6°—T 2 = (93) 



= S™ [4>, 5cf\ + Sh A f£S<t> a + ^Sh A 5h A + iSxAf^Tg [$\ 5u B 

+ i XAf*T 1 % fJ 5^6u J B + iSxa^T^S^uj 13 (94) 
The cubic and quartic terms are 

5(8+) = S (3+) + l5xA fA T a B ^ 6uJ B (95) 

Observe that r 2 is independent of the gauge fixing parameter, and that there are no h field lines. To take 
advantage of this fact, it is convenient not to couple sources to the h field. In this way, T 2 is independent of 
the h field, and the respective variations are exact, namely 

&AB~i[Gl 1 ] hhAB = 
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h x B 



We may then write the equations 



[ G X ] h4>Ap G 4>x + [G X ] hhAB Glx - $xhcS% 



as 



f^ x +^Gt r x = i5 xh c5 



KAC 



namely 



and 



riAr _ fA/~i0r 



X 



u hhA — J 



-fApG 



f3C 
<ph 



c 



-^fA^G^ + i5[ 



(96) 
(97) 
(98) 
(99) 

(100) 



which is of course what we expect from the N-L held being Gaussian. We could use these formulae to actually 
elliminate the N-L field from the 2PIEA, thus obtaining a reduced effective action. 
We shall assume that all fields with non zero ghost number vanish, i. e., 



u = x = G u 



G 



G u 



G 



x<t> 



G 



X h 







(101) 



Since the effective action itself has zero ghost number, it cannot contain terms linear on any of the above, 
and therefore this condition is consistent with the equations of motion. 

Given these conditions, we have, besides the equations determining the h propagators, the further equa- 
tions 



B-y 
T h<p 



leading to 



[ G %tf,a/3 G 0fc + t G \ 



[G 1 ]^ Q/3 + jfBafp 



[G 1 ] l p < p afj + jfBafp 

The second one gives nothing new, and the first yields 



G 



<phaB ^hh 
1 



■PC 
4>h 



fBaS 



BC 



[G 4,^/3 - G^ 
so finally we get the equation for the gluon propagator 



1 



S c ,a/3 — -^fBafp — i 



-fBafp 



ST 2 



a/3 6G 



a/3 



The other nontrivial equation is 



ST 2 



(102) 



In deriving this equation we must consider G^ uA and G^ %A as independent quantities. 



(103) 



(104) 



(105) 



(106) 
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3.3 The abelian case 

In the abelian case the classical action is quadratic, the 2PIEA is exact (T2 = 0). Since we know that 
S, a fiT% = 0, it is natural to decompose G af3 = G%? + Cf, where f B G%? = S, a0 G^ = 0. Then 

S, a0 Gp - ±f B J B ,pG^ = i5l (107) 

Observe that 

P£ f i = fB,{,[fB,-yTr\- 1 TS (108) 
satisfies P| = Pj, and kills vectors orthogonal to the gauge direction. In particular, P^S^p = 0, so 

and 

S, a pGp = i(62-P2 a ) (110) 

This shows that the dispersion relations in the transverse part are determined by the zeroes of the classical 
action, and are therefore gauge independent, unless (5% — Pl a ) is pathological. 

3.4 Some explicit formulae 

As before, the classical action is given by Eq. I|69(l . the quadratic terms by Eq. I|94[l and the cubic and 
quartic terms by Eq. 1)95(1 . where 

54,} = ±£$ y [$\ SFS^Sp + ^S% s drS^5^S4, s (111) 
The two loops approximation for T% reads 



r 2 = ^Xsis^sfs^) 



^ a f3y6\ u V "V "V "V / 2PI 

~s% [4>] s%,, [4>] (srs^5^sr's^'s^' /2pi 

c( 4 ) f-iaP^-fS , 1 ni(3) [7] ni(3) rri r aa' f~tfi$' f~m' 
' O a0-yS Lr <l>(l> Lr </><!> + Y2 Q/3 T ct'/3'7' L^J *0 U ^ °></> 

~nJa 1 lB0Ja' 1 \B'0' Kj i,<b u ujyA u ujvA' 



We may now write the Schwinger-Dyson equations 

= S c , a + h A f£ + -S c ^ n [4>] Gfy + ^f^Tf Ba G B uA + -^S^gSg^y [4>] G^G^P GJ^ 



- Sc, a fi-i[G 1 ]^ a0 + S^ jg G^ + -S^[^sf g ] y ,[$\Gf (t> G'^ > 

A_nfArpl fA'rpf' (~,B' S~tB 

= -itf'Tg [4>] + i [Gl 1 ]^ + l f^ B[j ff,Tf Bip ,Gf;G B ' xA (113) 
So finally we get the equation for the gluon propagator 
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S c ,a/3 ~ ~^fBaf/3 

eArpf fA'rpY 



GZ 



■ fArpl f/i rpj s~tB f~i 

l J~l 1 lBaJ-f' 1 lB'fl u ux^ U ' 



, c(4) r yS I o(3) r Tj c(3) 



55' /~tll' 



Goo s~i' 



(114) 



The other nontrivial equation is Eq. pi3|) . 

Let us assume there is a solution with = 0. Then, since S, a Tg — but also S„ [0] = 0, we have the 
identities 

S^isTb + S, a Ti B/3 = S^ a pTg [0] = 



(115) 



4 The Zinn-Justin equation 

We wish now to derive the Zinn-Justin equation appropiate to T2- The key observation is that under a 
BRST transform within the path integral, only the source terms are really transformed. Therefore 



The sources are 



(n [x r ]) j r + ^e rs e s (n [x r x s ]) K rs = o 



leading to 



(Sl[X r ]) 



K, 



ST 
5T 



-29 s 



5X r 
S „ 



- K. rs X" 



5G r 



- -e rs 9 s {SI [X r X s ]) K r 



SX r 
ST 

JXr 



<n [x r ]) x s - - (n [x r x s ]) 



[{SI [X r X s ]) -2 {SI [X r ])X s ] 



K, 



SG r 



(116) 



(117) 



(118) 



For simplicity, we shall assume that all background fields vanish. 

Since the Z-J operator has ghost number 1, it makes no sense to assume that all quantities with non 
zero ghost number vanish, as we have done in the previous section. However, we may still "turn on" these 
quantities one by one, and thus obtain partial Z-J identities. For example, we get three identities relating 
quantities with zero ghost number by requiring that the coefficients of u> and G^ vanish (we shall not 
investigate the first, as we are assuming no non zero backgrounds, and we are working throughout with the 
reduced 2P1EA). This means that we may still set 



to — x — G^uj — G xx — G x — G x h — 



(119) 



and retain only terms linear in G u <t> and G u h- In this approximation, terms with ghost number neither or 
1 must vanish identically, so 



15 



{n [lu d ] ) = <n [u, D u; E ] ) = {n [h A u, D ] ) = (n [<^ c ] ) = (o \xaxb}) = o (120) 

and 

5F ^ ^ ^ ^ =0 (121) 



^ «JG£? SGZ A 5G™ SG xxAB 

Also, since there are no preferred directions in gauge space, objects with a single gauge index must vanish 
out of symmetry, and therefore 

<n [</><*}) = (n [h A ] ) = (o [ XA ]) = §- [o] = (122) 



Finally, observe that at zero external sources, 

ST 5T ST 5T 



8h A SG^ hB SGh X AB 5G hh AB 



(123) 



In other words, from the terms in Eq. (jll8[l we keep the terms in <fxj>, fix, and \lo only. 

Eq. I|118|) must vanish at the physical point, since each coefficient vanishes. What is remarkable is that it 
vanishes identically, even if G$^Q. Now 8T/5G$ and 5T/5G B xA have ghost number zero, and therefore 
contain no terms linear in Ci A . We conclude that, to linear order in G°t£, we may write 



sr 

Where w means up to terms proportional to the equations of motion 
The transforms involve cubic terms 



WX^TT^A-O ( 124 ) 



Sl[J> a XA]=T%[(l>]u> B XA + i<i> a hA (125) 

The corresponding expectation values may be computed by adding to the classical action a new BRST- 
invariant source term 

R A x Jl [^Xa] (126) 

so we can write explicitly the action 

S = S c [4>] + h A f A r + ^h A h A + iXAftn m " B + +~4 xa [T% [4>] u b X a + icj> a h A ] (127) 
The quadratic terms are 

= [0] + h A f A r + ^h A h A + iXAftn [0] u B + n A xa [Tg [0] lu b X a + i4> a h A ] (128) 
The cubic and quartic terms are 

S(3+) = <^+) [4>] + ixa^T? b ^uj b + R A xa T^u B X A (129) 
so the missing expectation value is (cfr. Eq. I|78|l) 

(fi [^Xa]) = G B xA T% [0] + iG% hA + -JJ- (130) 



and the Z-J identity reads 

ST 



ST 

Guj xA Tg [0] + iG% lA + 



X A^B L"J T ™<f,hA T c-A 



xr< aA 



(131) 



Now 
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_6T_ 

zr<uA 
dG 4>x 



(132) 



and 



G, 



G B J [G-] 



XuiAB 



g: 



a/3 



G %[ G r%„ac < 1:5:5 1 



Tg[0] 



g f i^R IxujACJP + l^R Ix^C X R A 



a/3 



G C J « 



(134) 



We must still compute the derivatives w.r.t the external sources. Within the two loops approximation, 
we have 



{{xfS^ ds ^ d )T- Ep ^ e X a) 



fFrpJ rpa s~iE S~lD /~t&P 

2PI ~ J-y 1 lDS 1 lEp lj x^F X^A u 4„p 



so finally 



T 10 Tnl fFrp~f rpa r^E r^op f^-i — ll fA 

1 C M - h 1 lCS 1 lEp { - r x ^F t ~ r < t >< t > ^00 a - I L^fl J X u>AC J0 



a/3 



G C J « 



(135) 



(136) 



5 Gauge dependence of the propagators 

To investigate the gauge dependence of the 2PIEA, recall Eqs. (|35|l . 136|) and 13 7|) . Upon a change SF in 
the gauge fermion F, holding the background fields constant, we get 



With the same argument as for the 1PIEA, this leads to 

1 



ail*™ 



iUSFQ [X r ]) J r + -j6 rs (-r {,VF<> [X r X-]} K, 



Use Eq. (|TT7|) to get 

<5T|xr G ™ = i ((JFfi [JT]) ^- + [(<SFfi [X r A s ]) - 2 (*Ffi [X r ]) X s ] [ 



6G T 



(137) 



(138) 



(139) 



As before, we shall assume that all background fields vanish and that at such a point T^ r vanishes 
identically, so the above expression simplifies to 



6T\ Xr>Grs = Y rs ^T; Y rs = i (SFSl [X r X s ]) 
At the physical point, the Schwinger - Dyson equations now read 



(140) 



-r + y t 



SG tu 5G tu 5G rs 
Of course, the solution is now G tu + SG tu , so 

S 2 



SG tu SG r 



r = o 



r [SG rs + Y rs ] = o 



(141) 



(142) 



Since the Hessian is suppossed to be invertible, we must have SG rs = —Y rs . 

Let us also assume that all propagators with non zero ghost number vanish. Thus we are only concerned 
with 
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ST\ X r iG „ = {SFn [r<f] > + 2 (SFQ [ XA u D ] ) 



SF = -ixa { Sf A [0] + -5£h' 



(143) 



As before, we compute the corresponding expectation values by adding suitable sources to the action. 
These sources correspond to four and five - legged vertices. Now, recall the conventional argument that 

I - 1 = i - ^2 v n 

2i-^2nv n = (144) 
In our case, we have vertices with three, four or five legs, and I — 2, so 

2i - 3v 3 - 4v 4 - 5v 5 = 
i = 1 + V3 + Vi + v 5 

2-v 3 - 2v 4 - 3v 5 = (145) 

We are interested in the case where v§ = or 1. In the latter, we get an impossibility, so we must have 
t>5 = 0. Since we may discard the case when V4 = too, we must have v± = 1, V3 = 0, i = 2. In other 
words, to two loops, the only contributions to the expectation values are those where there are only four 
fields involved, and these contract among themselves. 

In conclusion, 



SG a 4 = -TZ10] 

= -nn 

= -?a [0] 



We may also write 



+ {a^0) 



Sf B_&i f B 

•h 2£ 7 
UJ i 2£ 7 



SG d>laB ~ G A a S T A [0] 



5 fp ~ 2| J? 



G^b + (a ~ 0) 



(146) 



(147) 



This shows that the zeroes of the inverse propagator are gauge invariant. To lowest order, the Z-J 
equation 



r r a rni f* 1 r r a f 1 *^ r 1 



iSp 



G 1, 



a/3 £ 



gZ « 



(148) 



implies 
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SG <pla(3 ~ £ [ G r] X u,AC 



= 7 /fl«W +/B^/a 

This is the result we wanted to show. 



(149) 



6 Propagator structure 

We shall try and apply the foregoing to clarify the structure of the 2-loop propagators. Let us begin with 
the equation for the ghost propagator (cfr. Eq. (|113fl ) 



/, 



A' 



B 



Ja> 1 lB0 1 lB'8' U rh<p u ujxA 



This suggests defining 



po 



rpOL fjl rArjia rpa 'p r^B' S~iB f 

1 B \_<P\ — Ja' 1 lB' 1 1 lB'f3' t ~'<t><p U ujxA b uxCJ) 



which is a projection operator 



(150) 



(151) 



pa pP pa 

Now consider the Takahashi-Ward identity (cfr. Eq. H136H 1 



(152) 



rpa rri] rpa fArpa r~*B /^ilp 
1 B L U J 1 lB'/3'Ja' 1 lBj u x^A 4><t> 



aX £ 



[ G r] x^ab fx ~° 



(153) 



It becomes 



That is 



Therefore 



Lfi 



G 4>1_ aX ~ 1 [ G r] x^ab fxG^xcfp ~ 



pa 
^LB 



aX £ 



fxfAB « 



(154) 



(155) 



Gij, 



, + jfxfAa, P2a 



r -l 

T<f><f> 



7 A 



(156) 



Comparing with the gauge-dependence identity Eq. 1(149(1 we see that the transverse part 
gauge-fixing independent. Comparing with the equation of motion Eq. 1)114(1 we get 



^Trbtf, 



7 A 







S c ,aB i 



T<p<p 



p + S a%s G lM> + 2 S °<h M S f&'i' [0] g h> g '11 



A fArpJ fA'rpY siB' (~lB 



(157) 



The decomposition of the inverse propagators leads to a related decomposition of the propagators them- 
selves. We have 



f A f >s r x P — p 13 

J 1 JAX'~'00 — -f£ 7 



(158) 



Therefore 
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where 



and 



Observe that 



and so 



We now have 



1 B VP] — Ja> 1 lB~ / 1 lB'f3' { - r < p<p ^uixA 



fA\G TU = 



DO T a rC 



pa p/3 _ tu f C T f D . f C T _ srC 



r -l 

T<p<p 



+ -zf\fAa 
aX ^ 



(159) 
(160) 

(161) 
(162) 
(163) 

(164) 



T<t>4> 



aX 



A/3 
T4"t> 



St- P. 







a La 



7 Appendix: Grassmann calculus 



(165) 



To be definite, we shall collect several known results concerning Grassmann calculus, which are necessary 
for the evaluation of the 2PIEA. For more details, we refer the reader to the monographs by Berezin |21| . 
DeWitt H2 and Negele and Orland [23] • 

Let us consider a set of independent Grassmann variables We define the left derivative <9/<9£" from 
the properties {d/d£ u ,C} = 6%, dl/d£ u = 0. 

We also define Grassmann integrals 



dC C =6 UV , dC 1 = 



(166) 

Observe that if £ = arj, then d£ = a~ 1 dr) (this follows from J e?£ £ = J dr/ rf). This allows us to prove the 
basic Gaussian integration formula 



d n £, exp \ — fM ut f \ oc y/detM 



Let us now consider the more general expression 



-1 



where the #'s are themselves Grassmann. Then, since M must be antisymmetric, 

CM UV C - 2i0„f = (f 1 - i (M- 1 ) 1 ™ d w ) M uv (f - i {M- X ) vx X ) - 9 U (M" 1 ) 1 



(167) 



(168) 



(169) 



d n i exp \ —CM UV C + i6 v C \ cx VdetM exp I -9 U (M" 1 ) 9 V 



(170) 



By differentiation, we get 
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8=0 



J d n i CCexp | -±CM UV F | = j d n t exp { yfM„f + iO v ? 

oc VdrfM (M- 1 )™ (171) 
This also follows from an integration by parts 



J d n iCC expire M uv c} = (M" 1 )™ J d n a u M wx ^exp^CM uv C^ 



(172) 



Observe that integrating by parts a Grassmann variable does not change the sign. 

We are interested in quadratic forms involving both Grassmann and normal variables. Let X r — (x a , £ u ), 
where the x's are normal, and consider the expression :X r M. rs X s := x a H a pxP + £, U M UV £ V + 2£ v N va x a . We 
shall call the object 



-N„. 



N u0 M uv ) 

a supermatrix. Observe that H is normal and symmetric, M is normal antisymmetric and N is Grassmann. 
Suppose the supermatrix M has a (right) inverse M^ 1 with elements 



Then 



H°f N 
-N u ? M 



VOL 

uv 



(174) 



Therefore 



H aP H& + N va N vl 
H a f)N v P - N ua M 
N U0 H^ - M UV N V ^ 
N U0 N V ? + M UW M WV 



UV 



= 5 1 



8 V 



H af} = [H a0 + N av (M- 1 ) 



t TV T 

N, 



! {3u 



-1 



M uv = [M UV +N 0U (fT 1 )""^ j 
We wish to check that both representations of N vl3 are equivalent. This means that 

, - Fir. 



Multiply both sides by [H f 



lfS5 + Nf3 v (M- l ) vu N Su 
(H~ 1 ) (3a N au M uv 



Multiply both sides by M vz + N 0V (H~ 1 ) (3a N a 



Hp 5 + N Px (M- 1 )^ N Sy ] = (M" 1 )™ N 5 . 



(175) 

(176) 
(177) 
(178) 

(179) 
(180) 



21 



- (H~ 1 ) tia N az [Hps + Nfr (M-^Nsy] = (M" 1 )™ N Sw [m vx + N Pv {H^f" N a 
We must check that 



(H-'f a N az N tix (M-r V N Sy = (M-') wv N Sw N [jv (if" 1 )"" N, 

Which is true, because the TV's are Grassmann and M antisymmetric. 
By the way, the left inverse is given by 



M7 1 = 



Let us suppose we use a supermatrix to implement a change of variables 

X r = M r s Y s 

Spelled out in full, 

x a = H al3 y p - N va f] v 
C = N^y + M uv r] v 



From the first equation 



So 



It follows that 



y = H- 1 (x + N V ) 



£ - NH~ l x = (M + NH-i-N) t] 



d n xd m i = d n xd m (£- NH^x) = det (M + NH^N) 1 d n xd m r] 
= det (M + NH^Ny 1 det (H) d n yd m r] 

So, if we define 

sdetM = det [M Q/3 ] det 
then X = MY implies dX = (sdetM) dY. Observe that 

sdetM" 1 = (sdetM)" 1 

Indeed 

sdetM^ 1 = [det (if + NM^N)] ^ dot (M) 

so we need to show that 

det (M + NH^N)' 1 det (H) = [det (H + NM^N)] [det (M)] _1 



or 



det (l + M^NH-^N) 1 = det (l + H^NM^N) 



(i; 



(1' 

(1' 



(1' 
(!■ 
(!' 
(!■ 
(!■ 



This follows from the TV's being Grassmann. To show this, observe that we can allways diagonalizc 
and reduce M to 2 x 2 blocks, so we may assume that H^ 1 = h is an scalar, TV = (61,62) and 
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M- . ( \ •» ) (194) 

Observe that 

while NM~ 1 N = 2to0i0 2 , so everything reduces to 

(1 - m/10102) -2 = l + 2mh6» 1 6> 2 (196) 

which is true because only the linear term in the Taylor development of the left hand side survives. 
We can now prove the Gaussian formula 

J dX cxp j^XMxj oc (sdctM) _1/2 (197) 

Let us now compute 

, p , 4 fdXX r X s cxp\^-XMX\ , 
(X r X s ) = J - — ? r—L 198 

Observe that 

X r x s = X r (Ml 1 ) 3 * M tu X u = 9 rs e rt (M^ 1 )"* X r M tu X u (199) 

On the other hand 



d 
dX 1 



cxpj^XMxj = -± [M tu X" + 0*X u M ut ]exp|^XMxJ 



:: \| \ :<[ _!XMX| (200) 



Finally 



Puting all together 



This shows in particular that 



—X r F = SlF + rt X r — (201) 



(X r X s ) = 6 rs (M^ 1 ) = (M^ 1 ) (202) 



In (sdetM) = r rs (M^ 1 )" (203) 
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